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INTRODUCTION 

Detailed  design  issues  associated  with  liquid  rocket  engine  injectors  and  combustion  chamber 
operation  require  CFD  methodology  which  simulates  highly  three-dimensional,  turbulent,  vaporizing, 
and  combusting  flows.  The  primaiy  utihty  of  such  simulations  involves  predicting  multi-dimensional 
effects  caused  by  specific  injector  configurations.  SECA,  Inc.  and  Engineering  Sciences,  Inc.  have 
been  developing  appropriate  computational  methodology  for  NASA/MSFC  for  the  past  decade.  CFD 
tools  and  computers  have  improved  dramatically  during  this  time  period;  however,  the  physical 
submodels  used  in  these  analyses  must  still  remain  relatively  simple  in  order  to  produce  useful  results. 
Simulations  of  clustered  coaxial  and  impinger  injector  elements  for  hydrogen  and  hydrocarbon  fuels, 
which  account  for  real  fluid  properties,  is  the  immediate  goal  of  this  research.  The  spray  combustion 
codes  are  based  on  the  FDNS  CFD  code’  and  are  structured  to  represent  homogeneous  and 
heterogeneous  spray  combustion.  The  homogeneous  spray  model  treats  the  flow  as  a continuum  of 
multi-phase,  multicomponent  fluids  which  move  without  thermal  or  velocity  lags  between  the  phases. 
Two  heterogeneous  models  were  developed:  (1)  a volume-of-fluid  (VOF)  model  which  represents 
the  liquid  core  of  coaxial  or  impinger  jets  and  their  atomization  and  vaporization,  and  (2)  a Blob 
model  which  represents  the  injected  streams  as  a cloud  of  droplets  the  size  of  the  injector  orifice 
which  subsequently  exhibit  particle  interaction,  vaporization,  and  combustion.  All  of  these  spray 
models  are  computationally  intensive,  but  this  is  unavoidable  to  accurately  account  for  the  complex 
physics  and  combustion  which  is  to  be  predicted.  Work  is  currently  in  progress  to  parallelize  these 
codes  to  improve  their  computational  efficiency. 

These  spray  combustion  codes  were  used  to  simulate  the  three  test  cases  which  are  the 
subject  of  the  2nd  International  Workshop  on  Rocket  Combustion  Modeling.  Such  test  cases  are 
considered  by  these  investigators  to  be  very  valuable  for  code  validation  because  combustion  kinetics, 
turbulence  models  and  atomization  models  based  on  low  pressure  experiments  of  hydrogen  air 
combustion  do  not  adequately  verify  analytical  or  CFD  submodels  which  are  necessary  to  simulate 
rocket  engine  combustion. 

We  wish  to  emphasize  that  the  simulations  which  we  prepared  for  this  meeting  are  meant  to 
test  the  accuracy  of  the  approximations  used  in  our  general  purpose  spray  combustion  models,  rather 
than  represent  a definitive  analysis  of  each  of  the  experiments  which  were  conducted.  Our  goal  is  to 
accurately  predict  local  temperatures  and  mixture  ratios  in  rocket  engines;  hence  predicting  individual 
experiments  is  used  only  for  code  validation.  To  replace  the  conventional  JANNAF  standard 
axisymmetric  finite-rate  (TDK)  computer  code^  for  performance  prediction  with  CFD  cases,  such 
codes  must  posses  two  features.  Firstly,  they  must  be  as  easy  to  use  and  of  comparable  run  times  for 
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conventional  performance  predictions.  Secondly,  they  must  provide  more  detailed  predictions  of  the 
flowfields  near  the  injector  face.  Specifically,  they  must  accurately  predict  the  convective  mixing  of 
injected  liquid  propellants  in  terms  of  the  injector  element  configurations. 

HOMOGENEOUS  SPRAY  COMBUSTION  MODEL 


The  homogeneous  spray  combustion  CFD  codes  utilize  very  general  thermodynamics  in  a 
conventional  CFD  code.  The  heterogeneous  codes  (described  in  Part  2)  use  tabulated  properties  for 
the  liquid  phase  and  ideal  gas  properties  for  the  vapor  phase.  Thermal  and  caloric  equations  of  state, 
vapor  pressure,  heat  of  vaporization,  surface  tension,  and  transport  properties  are  modeled  with  the 
equations  of  state  proposed  by  Hirshfelder,  et  al^  (we  term  these  the  HBMS  equations  of  state)  and 
with  conventional  correlations,^  for  the  other  properties.  The  property  correlations  used  were  not 
chosen  for  their  absolute  accuracy,  but  for  their  validity  over  a wide  range  of  temperatures  and 
pressures  and  for  requiring  a minimum  of  data  to  describe  a particular  species.  These  correlations 
are  explicit  in  density  and  temperature. 


HBMS  thermal  equation  of  state: 
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These  equations  are  based  on  the  "theorem  of  corresponding  states"  for  real  fluids,  which  essentially 
means  that  the  p-v-T  relations  for  all  species  are  similar  if  these  variables  are  normalized  with  their 
values  at  the  critical  point,i.e.  if  reduced  values  are  used.  The  reduced  values  in  these  equations  are 
indicated  with  a subscript  r.  Ho  is  the  ideal  gas  species  enthalpy.  Zc  is  the  compressibility  for  a given 
species  at  the  critical  point.  The  HBMS  equations  are  attractive  to  use  because  arbitrary  correlations 
for  vapor  pressure,  heat  of  vaporization,  and  liquid  densities  can  be  used.  Since  multi-component 
fluid/vapor  mixtures  may  be  present  in  the  flowfield,  the  mixture  properties  are  calculated  by  the 
additive  volume  method.  This  means  that  multiphase  mixtures  are  treated  as  ideal  solutions.  For 
H2/O2  propellants  under  conditions  where  the  species  become  ideal  gases,  the  thermodynamic  data 
from  the  CEC  code^  were  used. 

The  combustion  reactions  used  in  the  simulations  reported  herein  are  shown  in  Table  1.  Not 
all  of  the  reactions  were  used  in  all  of  the  combustion  simulations.  Elementary  rate  data  for  these 
reactions  are  reported  by  Gardner,  et  aP'*.  Such  data  are  empirical  and  were  obtained  for 
hydrogen/air  combustion,  under  conditions  far  different  from  those  encountered  in  rocket  engines. 
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Table  1.  Combustion  Model  for  H2/O2  Reaction 


Chain  initiation: 

H2  + O2  = 20H 

1.86  H2  + O2  = 1.645  H2O  + 0.067  O + 0.142  H + 0.288  OH 
Chain  Branching: 

H2  + OH  = H2O  + H 
2 OH  = H2O  + O 
H2  + O = H + OH 

O2  + H - O + OH 

Chain  termination: 

0 + H + M = 0H  + M 
20  + M = 02  + M 
2H  + M = H2  + M 

OH  + H + M - H2O  + M 


The  CFD  solver  used  was  the  Finite-Difference  Navier-Stokes  code  with  provision  for  using 
real  fluid  properties,  the  FDNS-RFV  code.  This  code  is  pressure  based;  it  differs  from  an  ideal  gas 
code  in  the  methodology  used  to  relate  the  pressure  correction  to  the  continuity  equation  and  of 
course  in  the  properties  subroutines  used.  The  pressure  correction  (p^)  equation  used  in  the  FDNS- 
RFV  code  is: 


p')  = -A*(p  Ui)  - 

At 

; Ui~  -Dp  ^p' 

where  the  superscripts  * and  n denote  the  value  at  the  intermediate  and  previous  time  steps, 
respectively.  Dp  is  the  inverse  of  the  matrix  of  the  coefficients  of  the  convective  terms  in  the  finite- 
difference  form  of  the  inviscid  equations  of  motion.  This  is  not  an  obvious  definition,  but  is  one 
which  has  made  the  FDNS-RFV  code  a useful  solver.  The  sound  speed  used  in  the  pressure 
correction  equation  is  that  calculated  for  the  real  fluid  multi-component  mixture. 

In  all  cases  simulated,  a k-E  turbulence  model  was  used  to  close  the  mass  averaged  transport 
equations  solved  by  the  code.  Our  experience  is  that  this  incompressible  turbulence  model 
overestimates  the  mixing  in  a combusting  flowfield.  However,  since  the  liquid  propellants  are  also 
mixed  by  this  model,  we  concluded  that  there  are  currently  insufficient  data  to  better  tune  the 
turbulence  model.  The  homogeneous  spray  model  has  been  used  to  simulate:  (1)  a single  element 
like-on-like  (LOL)  impinger  injector  element  and  a single  element  unlike  impinger  element  for  the 
configuration  and  flow  conditions  used  in  the  cold-flow  experiments;  (2)  an  ensemble  of  injector 
elements  in  the  Fastrac  engine;  and  (3)  several  configurations  of  the  vortex  engine  currently  being 
developed.^ 


Pp  P 
AT 


^•(uiPpP'  ) - p DpA 


P =P  +P  . P, 


= 7/ 


3 


RESULTS 


The  super-critical  combustion  case,  RCM-3,  was  simulated  with  the  homogeneous  spray 
combustion  model.  Any  drops  present  will  be  highly  unstable;  therefore,  this  model  should  represent 
the  flow  rather  well.  Local  equilibrium  and  simplified  finite-rate  combustion  submodels  were  used 
and  the  results  for  the  two  simulations  compared  well.  More  detailed  combustion  submodels  were 
attempted,  but  proved  to  behave  too  poorly  for  successful  simulations. 

The  preponderance  of  super-critical  spray  combustion  models  which  have  been  reported  have 
been  extensions  of  sub-critical  models.  Such  models  encounter  a basic  problem  in  over  emphasizing 
the  role  of  surface  tension.  Since  surface  tension  is  zero  for  super-critical  conditions,  drops  should 
not  exist.  Although  such  drops  can  be  observed  experimentally,  they  are  extremely  unstable  and  do 
not  survive  very  long.  The  homogeneous  CFD  model  was  developed  to  account  for  the  major 
physical  effects  which  do  exist.  Namely,  the  large  density  and  momentum  differences  which  exist  in 
multi-phase  super-critical  flows.  Such  a model  allows  one  to  accurately  relate  the  inlet  conditions 
at  the  injector  face  to  boundary  conditions  for  the  CFD  simulation.  This  relationship  is  essential  to 
predicting  the  effects  of  injector  element  configuration  and  inlet  momentum  vector  on  the  convective 
mixing  and  cross  winds  which  occur  in  practical  rocket  engines.  Otherwise,  one  is  forced  to  use  the 
historical  method  of  creating  costly  experimental  data  bases  from  which  to  choose  designs. 

The  injector  configuration  and  flow  conditions  for  the  supercritical  combustion  of  the  RCM-3 
test  case  are  presented  in  Fig.  1.  This  is  uni-element  shear  coaxial  injector  with  LOX  and  GH2 
propellants.  The  numerical  simulation  was  conducted  with  some  simplification  because,  initially, 
detailed  information  was  unavailable;  such  as:  (1)  the  flare  of  LOX  injector  near  the  exit  was 
neglected;  (2)  the  injector  was  flush  at  the  chamber  head-end  instead  of  protruding  into  the  chamber 
because  the  outer  diameter  of  hydrogen  tube  and  distance  between  the  chamber  head-end  and  the 
injector  exit  were  not  known;  (3)  the  nozzle  was  not  included  because  of  insufficient  information 
about  the  chamber  tail-end  and  nozzle  geometry;  and  (4)  the  coolant  (later  found  to  be  helium)  for 
the  chamber  wall  was  not  included  because  its  flow  rate  and  properties  were  not  specified.  As  can 
been  seen,  the  chamber  pressure  (60  bar)  is  well  above  the  critical  pressure  of  oxygen;  hence,  the 
homogeneous  real-fluid  model  was  used  to  simulate  this  test  case.  A two-zone  mesh  system  (61x39 
and  301x101)  was  used  to  model  the  injector  section  and  the  combustion  chamber. 

The  combustion  reactions  in  this  high  pressure  experiment  are  expected  to  be  in  local 
thermodynamic  equilibrium  and  were  simulated  as  such.  To  demonstrate  the  methodology,  two 
finite-rate  simulations  were  also  made  with  a subset  of  the  reactions  in  Table  1.  The  single  global 
reaction  which  produces  radicals  as  well  as  water  provides  a good  estimate  of  the  temperature  field. 
Its  rate  was  set  to  attach  the  flame  near  the  injector  tip.  Since  the  radicals  are  not  rigorously 
simulated  with  the  single  reaction,  a second  finite-rate  simulation  was  made  with  the  2-body  reactions. 
Backward  reaction  rates  are  determined  with  equilibrium  constants.  A third  finite-rate/equilibrium 
model  was  also  tested.  The  finite-rate  effects  were  described  with  the  9 elementary  reactions  in  Table 
1 with  the  combustion  assumed  to  be  in  local  equilibrium  when  the  temperature  was  greater  than 
some  specified  value.  Temperatures  of  1000  and  1500K  were  used  for  this  switch  point.  Thus,  the 
finite-rate  effects  would  be  considered  near  the  injector  and  in  the  expansion  section  of  the  nozzle. 
This  combustion  model  did  not  require  an  extremely  tight  grid  near  the  injector  tip  to  hold  the  flame. 
The  simulation  of  the  RCM-3  case  with  this  model  was  very  similar  to  the  equilibrium  combustion 
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model  For  high  pressure  cases  such  combustion  modeling  is  essential  to  keep  the  computation  stable. 

The  chemistry  and  turbulence  models  used  in  our  simulations  do  not  make  use  of  probability 
density  functions  (PDFs)  because  most  of  the  shear  layers  formed  by  the  injector  element  should  be 
continuum.  The  only  regions  for  which  this  might  not  be  the  case  are  the  intermittent  edges  of  the 
shear  layers.  Pope'°  terms  these  regions  the  "viscous  superlayer".  The  thickness  of  these  layers  are 
inversely  proportional  to  the  Reynolds  number  to  the  0.75  power.  For  these  high  speed  coaxial  jets, 
they  should  be  very  thin. 

The  flow  predicted  at  the  injector  tip  is  shown  in  Figure  2.  The  radial  temperature  profiles 
predicted  at  several  axial  stations  are  shown  in  Figure  3.  The  axial  profiles  at  several  radial  locations 
are  shown  in  Figure  4.  The  temperature  and  oxygen  and  OH  concentration  profile  fields  are  shown 
in  Figure  5.  The  combustion  models  used  do  not  predict  chemiluminescent  OH,  which  might  be 
observed  in  the  experiments.  These  results  are  shown  for  the  equilibrium  combustion  model.  Results 
for  the  finite-rate  equilibrium  combustion  simulations  are  very  similar,  hence  they  are  not  shown. 
The  wall  temperature  distributions  for  all  four  cases  are  compared  in  Figure  6,  and  as  noted  the 
results  are  very  similar. 

CONCLUSIONS 

The  following  conclusions  were  drawn  fi"om  performing  CFD  simulations  of  the  three  RCM 
test  cases  for  the  2nd  IWRCM. 

1.  A homogeneous  and  a heterogeneous  spray  combustion  CFD  models  have  been  developed  to 
simulate  combustion  in  rocket  engines.  Since  neither  of  these  models  is  expected  to  be  accurate 
until  critical  parameters  are  evaluated  fi-om  test  data,  simulation  comparisons  to  the  MASCOTTE 
type  experiments  are  needed. 

2.  The  utility  of  either  CFD  model  cannot  be  determined  until  values  of  critical  parameters  are 
determined  and  efforts  to  optimize  the  computational  efficiency  of  the  models  are  performed. 

3.  Although  the  CFD  rocket  engine  models  provide  much  more  detailed  information  concerning  the 
vaporization,  mixing,  and  combustion  process,  their  place  in  the  design  process  is  yet  to  be 
identified.  Older  more  approximate  rocket  "performance"  models  are  difficult  to  displace. 
Furthermore,  every  physical  process  thought  to  be  present  in  the  engine  does  not  have  to  be 
modeled  to  create  a useful  design  code.  There  are  more  knobs  to  adjust  in  the  code  than  there 
are  experimental  data  to  justify  their  turning. 

4.  The  experiments  conducted  in  preparation  for  the  2nd  IWRCM  appear  to  be  a significant  first  step 
in  providing  test  data  valuable  to  CFD  modelers.  However,  blind  comparisons  of  CFD  model 
predictions  to  such  data  are  premature.  The  CFD  modelers  have  not  previously  had  sufficient  test 
data  properly  specify  the  many  assumptions  which  are  necessary  to  simulate  such  complex  flows. 

5.  Better  communication  between  analysts  and  experimenters  needs  to  be  accomplished.  Can  the 
modeler  simulate  the  experiments  which  are  being  performed?  Can  the  data  obtained  from  the 
experiment  critically  test  the  model? 
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Critical  Pressure 

5.04  MPa 

1.29  MPa 

Critical  Temperature 

1 54.6  °K 

33  °K 

Inlet  Temperature 

65  “K 

287  “K 

Mass  Flow  Rate 

0.1  kg/s 

0.07  kg/s 

Figiire  4.  Axial  Profiles  of  Mean  Temperature  at  Various  Radial  Locations  of  RCM-3. 
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Figure  5.  Temperature  and  Species  Concentrations  Near  the  Injector  of  RCM-3. 
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Figure  6.  Near  Wall  Temperature  Distributions  for  Various  Chemistry  Model  of  RCM-3. 
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